Magnetoresistance oscillations in two-dimensional electron systems under 
monochromatic and bichromatic radiations 
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The magnetoresistance oscillations in high-mobility two-dimensional electron systems induced by 
two radiation fields of frequencies 31 GHz and 47 GHz, are analyzed in a wide magnetic-field range 
down to 100 G, using the balance-equation approach to magnetotransport for high-carrier-density 
systems. The frequency mixing processes are shown to be important. The predicted peak positions, 
relative heights, radiation-intensity dependence and their relation with monochromatic resistivities 
are in good agreement with recent experimental finding [M. A. Zudov et al., Phys. Rev. Lett. 96, 
236804 (2006)]. 
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The fascinating phenomena of radiation-induced 
magnetoresistance oscillations (RIMOs) and zero- 
resistance states (ZRS) in high mobility two-dimensional 
(2D) electron systems fi'^'^'^ continue to attract 
much experimentalSi'^'aiiaaiii2a3ii4^a6^ and 
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attention. 

Under the influence of a monochromatic microwave 
radiation of frequency aj/27r, the low-temperature mag- 
netoresistance of a 2D electron gas (EG), exhibits pe- 
riodic oscillation as a function of the inverse magnetic 
field 1/B. The oscillation features the periodical appear- 
ance of peak- valley pairs around node points lu/lOc = j 
(j = 1, 2, 3, 4, • • •, LOc = eB/m is the cyclotron frequency). 
With increasing the radiation intensity the oscillation 
amplitude increases and the measured resistance around 
the minima of a few lowest-j' pairs can drop down to 
zero. The nature of this observed ZRS, which was con- 
sidered as the result of absolute negative resistance, 
remains an issue for experimental confirmation. A re- 
cent measurement of photoresistance under simultane- 
ous illumination of two radiation fields of different fre- 
quencies by Zudov et al}''- provided a plausible evidence 
to link the ZRS to the negative resistance. However, 
the empirical superposition relation between the bichro- 
matic and two monochromatic resistances, needs refine- 
ment and justificationr^^ The present letter reports such 
an examination with the balance-equation approach^ to 
magnetotransport for high-carrier-density systems. 

In this photon-assisted-scattering based theory^^liiSS 
the transport state of a 2DEG subject to a dc field 
Eo and a normally incident monochromatic radiation 
Ej = Eis sm{ijjt)+'Eiic cos{ujt), is described by the electron 
drift velocity vq + cos(wt) + Vg sin(a;i) and an electron 
temperature Tc- They are determined by the balance 
equations (8)-(ll) in Ref.H^. The linear (vq — > 0) mag- 
netoresistivity is expressed as 
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where q|| = {qx,qy), is the electron density, ?7(q||) is 



the impurity potential, ^ = ^(q|| ■ Vc)^ + (q|| • Vs)'^/uj, 
and n2(q||,r2) = T^n2(q[|,fi) stands for the differentia- 
tion of the imaginary part of the electron density corre- 
lation function with respect to its frequency variable. 

The n2(q||,ri) function of a 2D system in a magnetic 
field can be expressed as sums over all Landau levels: 



n2(qii,f^) = 2^E^'^''(^B9'/2)n2(?,?',f^), 
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where Ci,i+.m{Y) = + my^-^Y'^e-^ [LfjY)] ^ with 
L™(F) the associate Laguerre polynomial, Ib = \f\J\eB\ 
is the magnetic length, /(e) = {exp[(e — [i)lT^ + 1}^^ 
is the Fermi distribution function at electron tempera- 
ture Te. The density-of-states (DOS) of the Z-th Lan- 
dau level is modeled with a Gaussian form, ImG;(e) — 
-(V2^/r) exp[-2(e - SiY IT\ having a half-width F = 

1/2 

{SeoJca/nm^o) around the level center e; = liOc- Here 
TO is the carrier effective mass, no is the linear mobility 
at lattice temperature T in the absence of the magnetic 
field, and a is a semiempirical parameter. 

It can be seen from expressions (0) and ^ that, in 
the case of low electron temperature (To ^ ep, the 
Fermi level) and large Landau-level filling factor p, for 
any integer N satisfying |iV| ^ i/, n2(q||,ri -I- Nujc) = 
H2(q||,il). Thus, when iv ~ luJc [l = 1,2,....), the 
H2(q||,nw) function in ([T]) is independent of n, that 

E„'/n(On^(qii,«^) = n;,(q||,o),and 
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■J2ql\U{<H^)\'U',{<H^,0). (4) 



indicates that at cyclotron frequency and its har- 
monics, uj = IlOc{1 = 1,2,...), the difference of R%x{C) 
driven by monochromatic frequency uj fields of different 
strengths, can appear only through the difference of elec- 
tron heating contained in the n2 function.— Under the 
condition of To <C ep, the effect of electron heating on 
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FIG. 1: (Color online) The magnetoresistivity R%nc induced 
by monochromatic radiations of frequency uj/2tt = 31 GHz at 
lattice temperature T = 1 K, in a GaAs-based 2DEG with 
iVc = 2.4 X lO^'^ m-^ ^io = 2000 mVVs and a = 5. 



photoresistance is quite weak)21 therefore, at integer-Wc 
value of frequency oj, the magnetoresistivity with radi- 
ation is essentially the same as that without radiation, 
RxxiO — Rxx{0), or the photoresponse vanishes at cy- 
clotron resonance and its harmonics. 

This is clearly seen in Fig. 1, where we plot the cal- 
culated magnetoresistivity i?"^ of a GaAs-based 2DEG 
with = 2.4 X lOi^m-2, = 2000 mVVs and a = 5, 
irradiated by linearly polarized monochromatic waves of 
frequency lu /27r = 31 GHz having incident amplitudes 
Ei = 0,0.5,0.7,1.0,1.3 and 1.7V/cm at lattice temper- 
ature T = 1 K. The scattering is assumed due to short- 
range impurities and the material parameters used are 
the same as in Ref. 33. the curves, with and without 
radiation, cross at uj/uJc — 1,2,3,4,5,6 and 7. This fea- 
ture, being pointed out before theoretically, ^ ^'^"i^^i^^ has 
been confirmed by experimentsi^>^"^"i° It is quite general, 
independent of the polarization state of the radiation, the 
behavior of the elastic scattering, and other properties 
of 2DEG, and holds quite precisely up to rather strong 
radiation field even multiphoton processes play impor- 
tant role, as long as the the 2DEG remains in degenerate 
(Tc ^ ep). It thus provides a convenient and accurate 
method to determine the effective mass of 2DEG. 

Another feature which shows up when the radiation 
strength is modest that contributions from two and 
higher photon processes are negligible and RIMOs result 
mainly from single-photon-assisted processes: 

+j! (0 (qii ,u;) + Jl, io (qii , -^)] . (5) 

Since the DOS function decays rapidly when deviating 
from the center of each Landau level, the magnitude 
of n2(q||,w) reaches a minimum of almost zero when 
uj/lUc = + 5) = 1,2,3...). Therefore around these 



half-integer-Wc positions, the magnetoresistivity 

^-(^)^-]^E9'lt^(qii)l'^o(e)n^(qihO) (6) 

is only weakly dependent on radiation strength through 
Jq(^), i.e. half-integer-Wc positions are approximately 
node points for modest radiation. This can also been 
seen in Fig. 1 where almost all the resistivity curves cross 
at Lo/ujc — 3.5,4.5,5.5 and 6.5, and the curves of lower 
radiation strengths cross at ui/iVc = 1.5 and 2.5. 

Under normally illuminating bichromatic radiation 
Eii = EiisSin(wit) -t- Eiiccos(wit) and Ei2 = 
Eii2s sm{uj2t) + Ei2c cos(w2t), the transport state of a 
2DEG can be described by the electron drift veloc- 
ity oscillating at both base radiation frequencies, vq -I- 
vicCos(LJii) + vissin(a;it) + V2cCos{uj2t) + V2s sm{LU2t) , 
together with the electron temperature To. They are de- 
termined by the balance equations (5)-(10) in Ref.lssl. 
The linear magnetoresistivity is expressed as 
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where ^1 = ^(^11 ' "^i^)^ + (^11 ' ^is)^/^^i and ^2 = 
• V2c)^ + (q|| • V2s)2/w2. (0) indicates that bichro- 
matic R%x^'' contains, besides the weighted superposition 
of mono-wi and mono-Lt;2 contributions, ^0 (Ci)-RS(?2) + 
Jq{^'2)R%x{^i)i also mixing uji and uj2 contributions 
(term with ni = 712 = and terms with both ni, n2 ^0). 

In the degenerate and large-r^ case, at positions U02 = 
lujc {I = 1, 2, 3....), n2(q|| , niLOi -1-712^2) function does not 
depend on 712, thus ([7]) reduces to 

qy 

00 

X E ^',(6)ni(q||,«i'^i)«i?^i(ei), (8) 

n 1 — — 00 

i.e. bichromatic R%x"^ {£,1 , £,2) equals the monochromatic 
R%x{£,i) under wi-Eii radiation (effect of electron tem- 
perature difference on photoresistance is negligible). This 
means that all the bichromatic R%x'^'^ {S,It^2) curves with 
fixed Eji but changing Ei2, cross with the monochromatic 
i?^^(fi) curve of Eii at integer-tJc points of 1^2. Likewise, 
at integer-Wc points of wi, ~ R%l{^2)- 

In the case of L02 ^ 1.5aJi, even and odd lOc points ex- 
hibit somewhat different behavior. At an even-Wc point 
of LOi {uji = 2ujc, 4a;c or QuJc), ^2 is also close to an 
integer ujc (^2 ~ i^c, and Qujc) that R%l{i2) is al- 
most identical to the dark resistivity Rxx{^)i independent 
of Ei2. Therefore, all bichro-resistivity i?"^'^^, mono- 
resistivity R%1. and mono-resistivity i?"^, cross at these 
points with dark resistivity Rxx{Q)- At an odd-ujc point 
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FIG. 2: (Color online) The magnetoresistivity Rxx induced 
by bichromatic or monochromatic radiations (Rxi'^^ , Rxx or 
Rxx) in the same 2DEG as described in Fig. 1. 
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FIG. 3: (Color online) Same as figure 2 in an enlarged scale. 



of LUi (wi = LOc, Swc, SoJc or 7lUc), UJ2 is close to a half- 
integer LOc {oJ2 ~ 1.5wc, 4.5a;c, 7-hujc or lO.Swc). It is a 
node point of R%'i{£,2) when Ei2 changes modestly. Thus, 
at these points bichromatic Rxi"^^ and monochromatic 
R'^^ and are nearly the same, but may be somewhat 
different from the dark resistivity. 
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Figure 2 demonstrates the calculated magnetoresis- 
tivity R%1.^'^ versus the inverse magnetic field for the 
same GaAs-based 2D system as described in Fig. 1 un- 
der linearly polarized bichromatic radiation of ti>i/27r = 
31 GHz and 102 = 47 GHz with incident amplitudes 
E-a = 1.5V/cm and = 1.5,1.8,2.0, and 2.3V/cm at 
T = IK. The monochromatic at E^ = 1.5V/cm 
and R%1. at Ei2 = 2.3V/cm are also shown. Comparing 
with monochromatic ones, the bichromatic curve consists 
of much more sizable peak- valley pairs. Since only minor 
high-order structures showing up in both monochromatic 
curves at these radiation strengths, most of the pairs ap- 
pearing in the bichromatic curves relate to mixing photon 
processes. We have identified the lowest order (|ni| < 1 
and 1^2! < 1) individual {uji or 072) and mixing (cji — 102 
or u)i +IJJ2) photon- assisted electron transition processes 
to all the pairs in the figure. Except those pairs around 
= 1 and LOi/uJc = 1, which relate to single-tJ2 and 
single-wi processes, and those around wi/wc = 2,4 and 
6, which relate to both individual and mixing photon 
processes, all other pairs are related to mixing processes 
uji+ UJ2 '■ I (I = ^1 2, ...17) with electron transitions jump- 
ing I Landau-level spacingsi^ 

A lower magnetic-field portion of this figure is re- 
drawn in Fig. 3 with enlarged scale for clearer compar- 
ison with experiments. The predicted positions and rel- 
ative heights of bichromatic peaks pi , p2 and pa around 
LOilujc — 2.7,3.1 and 3.9, and the predicted approximate 
crossing points ci (at W2/wc = 7), C2 (at wi/wc = 5) 
and C3 (at UJ2/0JC = 8), well reproduce^^ the experi- 
mental bichromatic peaks around 270, 222 and 182 G 
and the experimental crossing points marked in Fig. 2 
of RcLTt', where the empirical superposition relation is 
considered valid. The portion of ujc/^i = 1 vicinity of 
Fig. 2 is also replotted in the inset of Fig. 3 as func- 
tions of i3. A bichromatic peak p4 shows up clearly 
around ojc/oJi = 1.03, where mono-wi curve has a large 
positive value and mono-W2 curve exhibits a small posi- 
tive peak. This is exactly what observed experimentally 
around B « 730 G, where the superposition relation is 
apparently invalid. The present theory, instead, is able 
to provide a unified description for the bichromatic resis- 
tivity and its relation to individual monochromatic ones 
over the wide magnetic field range. 

This work was supported by Projects of the National 
Science Foundation of China and the Shanghai Municipal 
Commission of Science and Technology. 



* R. G. Mani, J. H. Smet, K. von Klitzing, V. Narayana- 
murti, W. B. Johnson, and V. Umansky, Nature 420, 646 
(2002). 

5 M. A. Zudov, R. R. Du, L. N. Pfeiffer, and K. W. West, 

Phys. Rev. Lett. 90, 046807 (2003). 
^ C. L. Yang, M. A. Zudov, T. A. Knuuttila, R. R. Du, L. 



4 



N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 91, 096803 
(2003). 

S. I. Dorozhkin, JETP Lett. 77, 577 (2003). 
^ R. G. Mani, J. H. Smet, K. von Klitzing, V. Narayana- 

murti, W. B. Johnson, and V. Umansky, Phys. Rev. Lett. 

92, 146801 (2004). 
^ R. L. Willett, L. N. Pfeifler, and K. W. West, Phys. Rev. 

Lett. 93, 026804 (2004). 
^° R. G. Mani, Physica E (Amsterdam)25, 189 (2004); Appl. 

Phys. Lett. 85, 4962 (2004). 
" M. A. Zudov, Phys. Rev. B 69, 041304(R) (2004). 

S. A. Studenikin, M. Potenski, A. Sachrajda, M. Hilke, L. 

N. Pfeiffer, and K. W. West, IEEE Tran. on Nanotech. 4, 

124 (2005). 

S. L Dorozhkin, J. H. Smet, V. Umansky, and K. von Kl- 
itzing, Phys. Rev. B 71, 201306(R) (2005). 

" R. G. Mani, Phys. Rev. B, 72, 075327 (2005). 

^5 J. H. Smet, B. Gorshunov, C. Jiang, L. Pfeiffer, K. West, 
V. Umanksy, M. Dressel, R. Meisels, F. Kuchar, and K. 
von Klitzing, Phys. Rev. Lett. 95, 116804 (2005). 

1^ M. A. Zudov, R. R. Du, L. N. Pfeiffer, and K. W. West, 
Phys. Rev. B 73, 041303(R) (2006). 

^"^ M. A. Zudov, R. R. Du, L. N. Pfeiffer, and K. W. West, 
Phys. Rev. Lett. 96, 236804 (2006). 

A. V. Andreev, L L. Aleiner, and A. J. Millis, Phys. Rev. 
Lett. 91, 056803 (2003). 

A. C. Durst, S. Sachdev, N. Read, and S. M. Girvin, Phys. 
Rev. Lett. 91, 086803 (2003). 
^° J. Shi and X. C. Xie, Phys. Rev. Lett. 91, 086801 (2003). 
X. L. Lei and S. Y. Liu, Phys. Rev. Lett. 91, 226805 (2003); 



X. L. Lei, J. Phys.: Condens. Matter 16, 4045 (2004). 
V. Ryzhii and R. Suris, J. Phys.: Condens. Matter 15, 
6855 (2003). 

2^ M. G. Vavilov and L L. Aleiner, Phys. Rev. B 69, 035303 

(2004) . 

'^'^ S. A. Mikhailov, Phys. Rev. B 70, 165311 (2004). 

J. Dietel, L. L Glazman, F. W. J. Hekking, and F. von 

Oppen, Phys. Rev. B 71, 045329 (2005). 

M. Torres and A. Kunold, Phys. Rev. B 71, 115313 (2005). 

I. A. Dmitriev, M. G. Vavilov, L L. Aleiner, A. D. Mirlin, 

and D. G. Polyakov, Phys. Rev. B 71, 115316 (2005). 

J. Earrea and G. Platero, Phys. Rev. Lett. 94, 016806 

(2005) ; Phys. Rev. B 72, 193414 (2005). 

28 X. L. Lei and S. Y. Liu, Phys. Rev. B 72, 075345 (2005); 

Appl. Phys. Lett. 86, 262101 (2005); 88, 212109 (2006). 
^° V. Ryzhii, Jpn. J. Appl. Phys., Part 1 44, 6600 (2005). 

C. Joas, J. Dietel, and F. von Oppen, Phys. Rev. B 72, 

165323 (2005). 

T. K. Ng and Lixin Dai, Phys. Rev. B 72, 235333 (2005). 
3^ X. L. Lei, Phys. Rev. B 73, 235322 (2006). 
^'^ A. Kunold and M. Torres, cond-mat/OeOSeW, 
^•■^ X. L. Lei and C. S. Ting, Phys. Rev. B 32, 1112 (1985); 
X. L. Lei, J. Appl. Phys. 84, 1396 (1998). 
This conclusion is in contrast with that of high tempera- 
ture and low Landau-level filling case. See X. L. Lei and 
S.Y. Liu, Appl. Phys. Lett. 82, 3904 (2003). 
If LOc/i^i = 1.03 is considered at 730 G, the pi, p2 and ps 
locate respectively at 263, 229 and 182 G. 



